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Abstract 

We argue that charmonium decays into meson pairs fall into two distinct 
classes: one that is under control of perturbative QCD and another one that 
is governed by a soft mechanism. We concentrate on a systematic analysis 
of J/ip(^f') decays into a light pseudoscalar and a light vector meson and rj c 
decays into a pair of light vector mesons. These processes belong to the second 
class and are characterized by non-conserved hadronic helicity. It is assumed 
that, in these cases, the charmonium state decays dominantly through a light- 
quark Fock component by a soft mechanism which is characteristic of OZI-rule 
allowed strong decays. Estimating the light-quark admixture by meson mixing, 
we obtain a reasonable description of the branching ratios for these processes. 
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1 Introduction 



Exclusive charmonium decays have been investigated within perturbative QCD by many 
authors, e.g. BL 0, H. It has been argued that the dominant dynamical mechanism is cc 
annihilation into the minimum number of gluons allowed by color conservation and charge 
conjugation, and subsequent creation of light quark-antiquark pairs forming the final state 
hadrons. This factorization of long- and short-distance physics has been shown to hold in 
the formal limit m c — ► oo. The dominance of annihilation through gluons is most strikingly 
reflected in the narrow widths of charmonium decays into hadronic channels in a mass region 
where strong decays typically have widths of hundreds of MeV Q . Since the c and the c quarks 
only annihilate if their mutual distance is less than about l/m c (where m c is the c-quark 
mass) which is smaller than the non-perturbative charmonium radius, and since the average 
virtuality of the gluons is of the order of 1 — 2 GeV 2 one may indeed expect perturbative 
QCD to be at work although corrections are presumably substantial. The charm-quark mass 
is too small in order to suppress power corrections decisively although it is large enough 
to allow perturbative QCD (pQCD) calculations. The bottomonium system, for which no 
exclusive hadronic decay has been observed as yet, should exhibit the pattern of perturbative 
predictions much cleaner. 

In hard exclusive reactions higher Fock state contributions are usually suppressed by 
inverse powers of the hard scale, Q, appearing in the process {Q = m c for exclusive char- 
monium decays), as compared to the valence Fock state contribution. Hence, higher Fock 
state contributions are expected to be negligible in most cases. For exclusive charmonium 
decays, however, the valence Fock state contributions are often suppressed for one or the 
other reason. In such a case higher Fock state contributions or other peculiar contributions 
such as power corrections or small components of the hadronic wave functions may become 
important. One such exception are the exclusive decays of the Xcj mesons for which the cc 
pair forms a color-singlet P-wave state in the valence Fock state (notation: cci( 3 Pj)). As 
has been shown in Ref. [|| the next-higher Fock state, ccg, where the quark-antiquark pair 
forms a ccs( 3 Si) state, is not suppressed as compared to the valence contribution. Therefore, 
the neglect of the ccg contributions, which are customarily referred to as the color-octet con- 
tributions ||, is unjustified in the Xcj decays. Indeed as has been discussed in Ref. |J for 
the PP channels (where P denotes a pseudoscalar meson) and in Ref. @ for their decays 
into baryon-antibaryon pairs the color-octet contributions are substantial and are definitely 
needed for obtaining agreement with experiment. 

Decays into final states involving vector mesons (V) have often peculiar properties too. 
On general grounds, only hadronic helicity non-conserving amplitudes contribute to some 
of these processes. Hence, these reactions are not under control of leading-twist pQCD. A 
famous example of the helicity non-conserving processes is set by the J /ip and VP' decays into 
pir. While J/tp — > pir has the largest branching ratio of all two-body hadronic J/ip decays, 
the branching ratio of — > pir is very small ||. The combination of these two experimental 
facts forms the so-called pir puzzle. Many of the other PV channels behave similarly. There 
are many attempts to understand the J/ip(*$>') — > PV decays; the proposed mechanisms reach 
from intrinsic charm in the light mesons Q, to color-octet contributions fpTof] , vector meson 
mixing JO, [l^] , final state interactions [ 13], 14 1 and admixtures of a glueball nearly degenerate 
with the J/ip (l5|, [l(|. Another example of a helicity non-conserving decay is set by the rj c 
decay into proton and antiproton. In Ref. [IT] it is attempted to explain this process through 
diquarks as quasi-elementary constituents of baryons. 
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In this article we are going to attempt a systematic study of the hadronic helicity non- 
conserving J/ip{^>') — ► PV decays. We will assume that, with a small probability, the 
charmonium possesses Fock components built from light quarks only. Through these Fock 
components the charmonium state decays by a soft mechanism. We will model this decay 
mechanism by J / ip — to — <p mixing and subsequent u (or <p) decay into the PV state. This, 
in the absence of the leading-twist perturbative QCD contribution, dominant mechanism is 
to be supplemented by the electromagnetic decay contribution and, depending on quantum 
numbers, by an anomalous, doubly Okubo-Zweig-Iizuka- (OZI-) rule violating contribution. 
We are not able to calculate the amplitudes of these decay processes directly but, employing 
flavor symmetry and mixing schemes for vector and pseudoscalar mesons, we find relations 
among them. With a few parameters, adjusted to experiment, we thus obtain a scheme of 
classification which comprises an understanding of these decay processes. The analysis of the 
J/ijj(^f') — > PV processes along these lines is partly to be considered as an update of previous 
work [l6|, 18, 19] where a general parametrization of the amplitudes for J/ip(^') — > PV 



is fitted to experiment. We improve this parametrization by including new ideas on meson 
mixing [2C] and by considering recent experimental results from the BES collaboration [21]. 
We also present a, partially even quantitative, physical interpretation of our parametrization 
which differs from those presented elsewhere [10, 1^, |l8|. Also new in our work is the extension 
of the mixing approach to the r] c — > VV decays. 

The paper is organized as follows: Qualitative features of charmonium decays into pairs of 
mesons are discussed in Sec. Next, in Sec. ^, we investigate the J/ip and decays into the 
PV channels by applying meson mixing. The mixing mechanism is also used to analyze the 
7] c — > VV decays (Sec. ||). In Sec. |B| we will comment on the mixing contributions to S-wave 
charmonium decays into baryons and antibaryons briefly before we present our concluding 
remarks in Sec. @. Details of our treatment of vector meson mixing are given in the Appendix. 



2 Qualitative features of charmonium decays into meson pairs 

Let us consider the PP, PV and VV final states. A few of these channels are strictly 
forbidden by angular momentum and parity conservation, see Table |l| Several other channels, 
characterized by non-conserved naturalness,[| a c / o"iO"2, are forbidden in pQCD to leading- 
twist order; i.e. higher- twist or other dynamical mechanisms are at work here. This comes 
about for the following reasons: The helicity amplitudes for these processes are proportional to 
the totally antisymmetric e-tensor contracted, in all possible ways, with the available Lorentz 
vectors, namely the two independent momenta, pi and p2, the polarization vector(s) of the 
light vector meson(s), £j, and, with the exception of the r] c , the polarization vector or tensor^Jof 
the charmonium state. Now, in the rest frame of the decaying meson, the polarization vector 
of a helicity zero vector (or axial vector) meson can be expressed as a linear combination of 
the two final state momenta regardless whether or not the mass of the light vector meson is 
neglected. Hence, the number of independent Lorentz vectors is insufficient to contract the 
e-tensor with the consequence of vanishing amplitudes for processes involving longitudinally 
polarized vector mesons. Moreover, ej(±l) = — £2(^1) holds in the rest frame of the decaying 
meson. As a consequence of these properties the only non-zero helicity amplitudes for the 
processes with non-conserved naturalness are those where Ai + A2 7^ 0. In other words vector 

1 The naturalness of a meson i having spin Ji and parity Pi is defined as at = (— l) Ji Pi 

2 Since the polarization tensor is symmetric only one Lorentz index can be contracted with the e-tensor. 
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Table 1: Charmonium decays into PP, PV and VV meson pairs. The symbols (— , e, y/) 
denote by angular momentum and parity conservation forbidden, in pQCD to leading-twist 
order forbidden and allowed channels, respectively. The brackets indicate that these channels 
violate either G-parity or isospin invariance for non-strange mesons. 

mesons are transversely polarized in reactions with non-conserved naturalness .[] We thus 
conclude that in these processes hadronic helicity conservation 

Ai + A 2 = (1) 

is violated. Helicity conservation is not a consequence of a particular symmetry but is a dy- 
namical consequence of leading-twist perturbative QCD (i.e. using leading-twist wave func- 
tions and valence Fock states only): The virtual gluons from the annihilation of the cc pair 
create the light, (almost) massless quarks and antiquarks in opposite helicity states. To the 
extent that the hadronic wave functions do not embody any non-zero orbital angular momen- 
tum components, the quark helicities sum up to their parent hadron's helicity. Hence, the 
total helicity of the final state hadrons is zero. Consequently, processes that violate helicity 
conservation, are not governed by leading-twist pQCD. The remaining two-meson channels, 
marked by ticks in Table [j], are accessible to a perturbative treatment. 

Next let us consider G-parity and isospin. G-parity or isospin-violating decays are not 
strictly forbidden since they can proceed through electromagnetic cc annihilation and may 
receive contributions from the isospin-violating part of QCD. The latter contributions, being 
of the order of quark mass differences, seem to be small Q. G-parity or isospin-violating 
decays of G-even charmonia (e.g. r/ c , Xcj — ► PV (J = 1, 2) for non-strange final state mesons) 
have not been observed experimentally as yet ||]. Proceeding on the assumption that these 
decays are dominantly mediated by cc — * 27* —* PV, this is understandable. They are then 
suppressed by a factor (a em /a s ) 4 as compared to the G-parity and isospin allowed decays 
of the G-even charmonia and their decay widths are therefore extremely small. Channels 
involving strange mesons (e.g. KK*), should also be strongly suppressed by virtue of [/-spin 
invariance. For J/tj) decays the situation is different. Many G-parity violating (e.g. ir + ir~) 
or isospin-violating (e.g. lott) decays have been observed, the experimental branching ratios 
being of the order of 10 -4 — 10~ 3 Q. As compared to G-parity and isospin allowed J/i/j 
decays they are typically suppressed by factors of about 10~ 2 — 10 _1 in accordance with what 
is expected for an electromagnetic decay mechanism^] (see Fig. p]). 

That at least two distinct dynamical mechanisms are at work in exclusive charmonium 
decays can be realized from a comparison of J /if) and \&' decays. Suppose these decays 

3 Somewhat exceptional is the Xci — * VV case where one out of the three spin-1 mesons may be longitudi- 
nally polarized. 

4 A possible correction due to cc annihilation mediated by 7*3*3* is ignored by us here. 
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Figure 1: Electromagnetic J ftp and ^f' decays into meson pairs. The shaded circle indicates 
the time-like electromagnetic (transition) form factor of light vector and/or pseudoscalar 
mesons. 

are under control of a mechanism that respects QCD factorization and where the charm 
and anticharm quarks annihilate into gluons and/or photons at small mutual distance. The 
charmonium wave function is then probed at small spatial separations and, therefore, is well 
represented by the corresponding decay constant, /ju or f^>, which is measured in electronic 
decays of that charmonium state. If this short distance mechanism is responsible for the J ftp 
and W decay into a particular two hadron channel, the ratio 

Bj&^hlM ffW^e+e-) ( Q 12 (J/ip) \ n ( ) 

Kl2 B(J/i> -» hM B(W -> e+e") I ei2 (*') J 1 ' 

should be close to unity^ (n = 3 for PV channels and n = 1 otherwise). Since B{^' — > 
e + e~)/B(J/tp — > e + e~) amounts to about 14%, the relation (g) is occasionally termed the 
14% rule. The phase space factor in Eq. (0) is defined by 

f?l 2 (j) = s[l ~ 2(mf + m|)/M? + (mf - m^/M/ (3) 

where Mj is the mass of the charmonium state. It is of numerical relevance only for particles 
with masses above 1 GeV. Experimental results for k are listed in Table 0; the data are 



taken from Refs. |8|, 21, 24, g^, |26|, 27]. For comparison we also include experimental results 
for the baryon-antibaryon, vector-tensor and axial-vector-pseudoscalar channels in the table. 
We observe that, within often large errors or only within bounds, k is indeed compatible 
with unity for the baryon-antibaryon channels, for b\ir, K*°K*° and for some of the PV 
channels, notably for the G-parity and isospin-violating channels ttuj, rjp, rfp. For most of 
the other PV channels and, perhaps, to a lesser extent also for the vector-tensor channels, 
K is well below unity. We stress that the decays into a vector and a tensor meson are 
not forbidden by hadronic helicity conservation, but leading-twist pQCD feeds the helicity 
amplitude Moox c only. There are many other, in general non-zero amplitudes that do not 
respect hadronic helicity conservation. Note that the branching ratios of J/ij) decays into VT 
are large, comparable with those for the PV channels while those for the corresponding VP' 

5 In a consistent perturbative calculation to lowest order in the c-quark velocity expansion j^, ^, ^] one has 
to evaluate the hard scattering amplitude from 2m c instead from the bound state mass; the mass difference is 
an effect of order v 2 , i.e. of the same order as the contributions from the next higher Fock states. An additional 
factor of (M t //Mj^) 2 " in Eq. (H) (n = 4 for baryon-antibaryon channels for instance), to be found in the 
literature occasionally, is inconsistent in this respect. 
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Figure 2: Schematic picture of the mixing mechanism for charmonium decays into meson 
pairs. 



decays are small. For the PP channels the situation is unclear; better data are demanded. 
We do not include the data for multi-particle channels in the table but we note that, with 
the exception of the final states consisting of three pseudoscalars, k is compatible with unity 
for all these channels. This fact is presumably a consequence of an underlying perturbatively 
generated three gluon jet structure which fragments into multi-particle final states. This 
interpretation is supported by estimates of the total width for J/ip decays into light hadrons 
through the decay into three real gluons [28]. In the case of the T decays the three jet 
structure is experimentally established [p9j] . 

The small value of k for those PV channels that are dominated by strong interactions, 
arises from large J/ip branching ratios as an inspection of Table [2| reveals. The \P' decays 
into PV, on the other hand, seem to behave, at least in tendency, according to expectations 
based on pQCD; their branching ratios are smaller than those for the leading- twist allowed 
channels. Obviously, another dynamical mechanism is called for which is active for the J/ip 
decays but suppressed for the ones. For this mechanism we assume that the charmonium 
state possesses Fock components built from light quarks only. It can then decay through 
these Fock components by a soft mechanism that is characteristic of OZI-rule allowed strong 
decays (see Fig. ^). This mechanism is obviously more peripheral than the short distance cc 
annihilation; i.e. it probes the charmonium wave function at all quark-antiquark separations 
and feels therefore the difference between a IS and a 2S radial wave function. The node in 
the latter is supposed to lead to a strong suppression of this mechanism in the \P' decays.^] 
We model this mechanism by mixing of vector mesons, J/ip — uj — (p. For the mixing angle, 9 C , 
being of order 1/mj?, we take the value —0.1° which corresponds to a light-quark admixture to 
the charmonium state of order of 10 -3 , see the Appendix. The light vector mesons, to and </>, 
couple to the final state mesons through vertex functions, e.g. g^vp- I n terms of vector meson 
mixing the suppression of the \P' — > PV decays follows from substantially weaker couplings 
(arguments in favor of that assumption are given in Refs. 12, |3l|]) of the radially excited u>' 
and <p' mesons (assumed to mix with the to the ground state pseudoscalar and vector 
mesons and, perhaps, from a smaller mixing angle. The mixing mechanism may also be at 
work in the vector-tensor channels (with vertex functions gvVT) an d may generate the small 
value of k there, see Table g 

Conversely one may also consider the cc Fock components of u and <p mesons in the final 

6 This assumption bears a resemblance to the node effect discussed in Ref. |30|. 
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state; their probabilities are the same as the ones for light quark components in the J ftp 
meson, see, Eq. (^Tj) in the appendix. The initial cc pair from the J/ip may then feed that 
Fock component of the u or <p. However, the second final state meson is to be generated by 
an OZI-rule violating mechanism now, and is to be considered as a higher-order correction 
to the leading mixing mechanism we propose. For channels like urf one may also think of 
an excitation of the cc components of both the light mesons. Obviously, the probability of 
finding a cc pair in a light meson enters quadratically now, and this decay mechanism is 
suppressed by an additional power of \jrr? c . Hence, we ignore both the above possibilities. 
The initial cc pair may also excite a ccqq Fock component in the light (vector) meson, as 
suggested by Brodsky and Karliner |J. In that picture the qq pair from this Fock component 
forms the other final state meson. Of course, one expects the probability of this higher Fock 
state to be very small. In any case the approach of Brodsky and Karliner should be viewed 
as complementary to ours, since the mechanism that they propose cannot be represented by 
meson-mixing supplied with OZI-rule allowed decay vertices. 



We note that in Refs. |10L 16, Qq] the strong interaction mechanism for the hadronic 



helicity non-conserving decays although parametrized in a similar way as we do, is interpreted 
differently. Thus, Chen and Braaten argue that this contribution is generated by an 
additional Fock state gluon (the cc pair is therefore in a color-octet state) which carries a 
unit of helicity from the J ftp to one of the light mesons without participating in the hard 
process. Chen and Braaten provide arguments that the suppression of the W decays into 
PV mesons as well as into the vector-tensor channels is due to the energy gap between the 
mass of the ^f' and the DD threshold. In Refs. |i~6| , |is} |, on the other hand, mixing of the 
J/ip with a fairly narrow J PC = 1 glueball, nearly degenerate with the J /if), is assumed. 
The J ftp can then decay through the glueball while this mechanism is not available for the 
5 r '. Searches for such a glueball, performed by the BES collaboration 32], however turned 
out negative although not fully conclusive to claim the demise of the glueball interpretation. 
For the final state interaction mechanism advocated in Refs. [13, 14] the suppression of the 
vp' — > PV decays is due to a fortuitous cancellation of various contributions and, probably, 
does not hold for the vector-tensor channels. 

The electromagnetic decay mechanism, proceeding through the subprocess cc — > 7* — > 
hili2 (see Fig. |l]), also probes the charmonium wave function at small spatial separation and, 
hence, would lead to k ~ 1 if it dominates. This is what we observe, or at least what is 
indicated, by the isospin-violating ttlv, rjp and rj p channels. For the other J/ip(p&') — > PV 
decays a superposition of the mixing and the electromagnetic mechanism occurs. For the 
flavor-singlet channels one has to consider an anomalous (formally doubly OZI-rule violating) 
contribution |l8| , |i"9f| too, see Fig. ||. The mixing mechanism also applies to the rj c — > VV 
decays with vector meson mixing replaced by that for pseudoscalar mesons pC|]. 

The Xcj may also decay through light-quark Fock components. However, in terms of 
meson mixing, the model would require the coupling of appropriate ++ , 1 ++ , 2 ++ mesons 
to pairs of ground state mesons. These couplings are expected to be small and, therefore, 
the Xcj — ► PP, W {J = 0, 2) decays seem to be dominated by perturbative contributions 
while the hadronic helicity non-conserving decays Xcl ~^ W should be strongly suppressed. 
Indeed, the latter decays have not been observed experimentally as yet ||. As we already 



mentioned in the introduction the perturbative analysis |5|, 23] of the Xcj decays into the 
PP channels provides results in fair agreement with experiment [33|] . In this analysis, it 
is however important to include the color-octet contributions, i.e. the ccg Fock state of the 
Xcj- The contribution of the valence Fock state alone, evaluated from the asymptotic form 
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Table 2: Exclusive J/ 1(1 and ip' branching ratios and their ratios scaled by the corresponding 
electronic branching ratios and phase space corrected according to Eq. (0). Data are taken 

|2l 



from |8j in the J /if) cases and from the BES collaboration 24 
otherwise. Some of the BES data are preliminary; the BES data on 
quoted in Ref. pi. 
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Figure 3: Schematic picture of the anomalous (doubly OZI-rule violating) decay mechanism 
for J/ip and decays into PV meson pairs. 



of the pseudoscalar meson's light-cone wave function (as determined from the analysis of the 



P7 transition form factors [34, 35, 36|), is substantially smaller than experiment. Aspects 

The extension of the 



of the 77 — 7/ mixing in these decays have been discussed in Ref. 
perturbative analysis to the VV channels, for which first data exist [33], is straightforward 
but the required wave functions of the vector mesons are not known to a sufficient degree of 
accuracy at present. 



3 Mixing mechanism and the decays J/ip(ty') — ► PV 

Let us now turn to the J/ip and decays into a pseudoscalar and a vector meson. According 
to the discussion presented in Sec. || we take into account the three mechanisms shown in 
Figs, [l], || and ||. Each of these contributions will be parametrized by a reduced amplitude, 
common to all PV channels, multiplied by flavor factors, meson-mixing factors and flavor- 
symmetry corrections. Because of the lack of a compelling theory we have to treat the reduced 
amplitudes as free parameters to be adjusted to experiment. Despite this fact our ansatz which 



is structurally similar, although different in detail, to those discussed in Refs. [[LC], 16, [L^, [19| 
provides a systematic description of the J /ip and VP' decays into PV mesons as we will see 
shortly. 

To begin with we discuss the electromagnetic decay mechanism shown in Fig. |l|. Decom- 
posing the helicity amplitudes for a decay of an n 3 Si charmonium state into PV covariantly 
as 

A4oav,a(« 3 Si-VP) = ffiffg^ 1 ? e(pi,p 2 ,e* v (\ v ),e(\)) , (4) 

we find for the electromagnetic contribution to the invariant amplitude 

^(n 3 Si) = 47ra em e c e^/(n 3 S 1 )M(n 3 S 1 )F P y( S ) . (5) 

/(n 3 Si) is the decay constant of the n 3 Si charmonium state. Its numerical value is taken 
from Ref. @ (f(J/i>) = 405 MeV and /(*') = 282 MeV). F PV is the time-like P -> V 
transition form factor for transversely polarized vector mesons at s = M 2 (n 3 Si). We write 
this form factor as 

Fpv(s) = cf v F p oAs) , (6) 
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Table 3: The flavor factors appearing in the amplitudes for n 3 Si — ► decays. For definitions 



see Eqs. (g,|l^,|l2|) . The parameters for rj — rj' mixing, <pp and 0g, are defined in Ref. [20], 
those describing vector meson mixing, (j)y, 9 y , y, in the Appendix. 



and treat ^0^0(5) as a free parameter (cy^o = 1 by definition), dp^y depends on the electric 
charges of the PV system in units of the elementary charge eo- We also absorb rj — 7/ and 
u — <j) mixing factors as well as a flavor symmetry breaking correction factor, y em , into Cpy. 
The parameter y em takes into account the difference between the occurrence of strange and 
non-strange quarks in the form factors; its value will be estimated below. The mixing of the 
pseudoscalar mesons is described in the quark-flavor scheme advocated in Ref. PQ| . Vector 
meson mixing is treated analogously, the details are presented in the Appendix. Since the 
vector meson mixing angle, <j>v, is so small (see Eq. (j30|)) we put cos^y ~ 1 and sin^y ot 4>v- 
The flavor factors &py are compiled in Table ||. In Eq. (||) we also allow for a phase, 9 e , 
relative to the strong interaction contributions discussed below. The transition form factors 
are therefore to be understood as absolute values. The invariant function Apy in Eq. (Q) is 
normalized in such a way that the decay width reads 

r(n Sl ^PV) = — M(n3gi) \A PV \ . (7) 

The four isospin- violating channels now allow a simple extraction of the p°ir° form factor. 
From a fit to the experimental branching ratios of the isospin-violating decay channels we 
find 

F p °AM 2 j lf ) = 27-10" 3 GeV" 1 , 
F pH o(Ml,) = 19 • 1(T 3 GeV^ 1 . (8) 

The results of the fit are listed in Table Note that the values (|8|) of the form factors are 
subject to errors which we estimate to 7.5% and 21% for the J/tp and cases, respectively. 
The calculation of the time-like form factor, in particular in the vicinity of resonances, is a 
delicate issue and, in so far, reliable theoretical information on this form factor is lacking. 
Since the p meson is transversely polarized, the form factor is under control of higher- twist 
contributions and should therefore fall off as 1/s 2 asymptotically. Comparison of the two 
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PV 


10 4 • B tt (J/V) 


10 4 • S cxp (J/^) 


10 4 • B th (*') 


10 4 • B cxp {V) 


LU1T 


3.78 


4.20 ±0.60 


0.38 


0.38 ± 0.20 


PV 


2.04 


1.93 ±0.23 


0.21 


0.21 ±0.12 


pv' 


1.06 


1.05 ±0.18 


0.12 


< 0.3 


(j)7T 


0.013 


< 0.068 


0.001 




pir 


127 


128 ± 10 


0.13 


< 0.28 


0J7] 


14.7 


15.8 ± 1.6 


0.12 


< 0.33 


UJT]' 


1.86 


1.67 ±0.25 


0.63 


0.76 ± 0.48 


K*+K- + c.c. 


48 


50 ±4 


0.13 


< 0.30 


K*°K° + c.c. 


42 


42 ±4 


0.51 


0.81 ±0.29 




7.5 


6.5 ±0.7 


0.16 


0.35 ±0.20 


H 


2.6 


3.3 ±0.4 


0.46 


< 0.75 



Table 4: Comparison of predictions and experiment for the branching ratios (multiplied by 
10 4 ) for J/tp and ^' decays into a pseudoscalar and a vector meson. The results for the first 
two sets of decay channels are obtained from fits (see text). The third set of decay channels 
involving strange mesons are predictions based on the parameters determined in the fits to 
the first two sets of data. The J/tp data are the PDG averages ||, the ^' ones are preliminary 



BES data [21, ^4|, 27]. The values of the mixing parameters for vector mesons are quoted in 



the Appendix, those for pseudoscalar mesons are 4>p = 39.3° and 9s = —21.2° [pQ |. 

values in Eq. (||) rather reveals an approximate 1/s behavior in the considered region of 
s. A recent QCD sum rule analysis [ |39"[ of the space- like pir form factor yields values that 
are somewhat smaller than those quoted in Eq. (||). The enhancement of time- like form 
factors relative to space- like ones is known from the pion and nucleon cases; the disparities 



typically amount to about a factor of 3 in the charmonium region (see e.g. [17, 40 1 ) . In 
any case, measurements of the PV transition form factors in the charmonium region would 
be welcome. Such data would permit a critical examination of the interpretation of the 
electromagnetic contribution (|5|) (cf. also the discussion in Ref. (l(|]). In particular it might 
be checked whether or not the neglect of contributions from the isospin-violating part of QCD 
and from cc annihilation mediated by "f*g*g* is indeed justified. 

The 4>tt channel offers direct possibility of testing uj — 4> mixing. As inspection of Table || 
reveals the predicted J/ip — > (fttr branching ratio obtained from our value of the mixing angle 
(p0|), is well in agreement with the experimental bound. A measurement of the branching 
ratio for the 4>tt channel would allow a direct determination of 4>y through the ratio 

B(n 3 S 1 -+ <M) 2 (QuJjf^Y , . 

a relation that was proposed by Haber and Perrier long time ago. 

The invariant amplitude for the strong decay mechanism (see Fig. |2|) is parametrized as 



yso = c ^«r ix (n 3 si) , (io) 



where dpy depends on flavor symmetry breaking and meson mixing factors. As for the 
electromagnetic contribution we put c™ lx Q = 1. By virtue of the node effect (see the discussion 
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in Sec. (2j) this mechanism does not contribute to the \&' decays, i.e. we assume g mix (J/tp) 3> 
g mix (^') ~ for the reduced amplitude. 

The anomalous contribution, depicted in Fig. can only contribute to the flavor singlet 
meson channels. It is assumed to proceed via a short-distance annihilation of the initial cc 
pair and a subsequent creation of the rj and rf states controlled by the matrix element of the 
topological charge density (Q\jzGG\rj^) which incorporates the U(1)a anomaly. As shown 
in Refs. [20, 37] the ratio of both the gluon matrix elements is given by 



(0\%GG\r/) 



tan6> 8 , (11) 



where the angle 6s is related to the rj — rj' mixing angle <p>p via 9$ = (j)p — arctan[\/2/y], 
where y being the ratio of basic decay constants, / q and / s in the quark-flavor scheme p0| . 
The numerical values for the mixing parameters, determined in Ref. [pC[| , read <pp = 39.3°, 
6s = —21.2° and y = 0.81. In view of these considerations we parametrize the anomalous 
contribution as 

Ap^in^Si) = f{n 3 Si)c P n $ m g anom {n 3 S 1 ) , (12) 

and define c^? m = y/2. The flavor factors c P n y m for the other channels can again be found 
in Table ||. Because of the short distance cc annihilation mechanism g anom ) the reduced 
amplitude of the anomalous contribution, should only mildly depend on the charmonium 
mass; we therefore assume g anom ( J/tp) = g anom (\I/ / ) for simplicity. 

A fit to the J/tp and ^f' decay widths into pit, urj and urf provides the results shown in 
Table || for the following values of the parameters 

g m[x (J/ip) = 0.020 GeV , 
g anom = -0.014 , 

6 e = 78° . (13) 

The value of the relative angle between the electromagnetic and the strong contributions is 
very similar to that found by Bramon et al. l£j. Relative phases with values close to 90° 
between different decay amplitudes have also been observed in other approaches to J/tp (i/j 1 ) 



decays [14, 42, [43|. Concerning the relative size of contributions from the three individual 
decay mechanisms (electromagnetic, mixing, anomalous), the result strongly depends on the 
considered decay channel. For instance, in the case of J/tp — * pn the ratio of [vl 6111 !/!^ 111 "! 
is only 0.1 whereas for J/tp — ► K*°K° the same ratio is about 0.3. The anomalous con- 
tribution tends to interfere with the mixing contribution destructively and is in some cases 
very important. This is most drastically seen in the channel J/tp — > (prf where the ratio 
\A em \/\A mix + vl anom | is about 1.1. 

The decays into KK* and <prj(') involve strange quarks which necessitates the consideration 
of flavor symmetry breaking. In the case of the electromagnetic contribution one has to 
pay attention to the fact that the PV transition form factors are controlled by higher-twist 
contributions since the vector mesons are transversely polarized. Hence, the form factors 



are proportional to a mass scale which, as shown by Efremov and Teryaev [44|, is of the 
order of the hadron masses and not, as one may naively expect, set by current-quark masses. 
Modelling this mass scale by appropriate effective constituent-quark masses (fh), we estimate 
ycm ^ ms/rhq ~ 1.3. For (prj and cj>rf channels the factor should appear quadratically. For 
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the mixing contribution strange quarks come either from J /if) — <f> mixing, and are therefore 
oc cot 8y, or by a soft creation of a ss pair out of the vacuum. This creation process is 
suppressed by a factor that is less than unity, relative to the creation of a non-strange quark- 
antiquark pair. In accordance with the flavor symmetry breaking in vector meson mixing, 
see the Appendix, we identify this factor with y and use the value 0.8 for it. Note that in 
our mixing scheme \/2cot^ = —y, see Eq. (pop. Therefore, the pattern of SU(3)f breaking 
is similar for the mixing and the electromagnetic contributions. With these values for the 
flavor symmetry factors and the other already determined parameters at hand we are in the 
position to predict the branching ratios for the decay channels involving strange quarks. The 
results are also shown in Table ||. The quality of the predictions for the entire set of J / ip and 
ty' decays into a PV pair is good. For the J /if) decays the quality is similar to that obtained 
in Ref. although we have a smaller number of adjustable parameters. 

In our approach the ^f' decays into the PV pair electromagnetically with the exception 
of the flavor-singlet channels where the anomalous mechanism contributes as well. In so far 
our analysis differs from that advocated by Tuan |l6|] who extended the analysis in Ref. jn| 
to the case of the Tuan assumes the electromagnetic and the strong contributions to 
cancel approximately in the ^' — > pit channel. This assumption leads to substantially smaller 
branching ratios for the W — > pit and \&' — > K* + K~ channels than we obtain. Since for 
both channels the BES collaboration only provides bounds the issue which of the contending 
models is correct, cannot be settled as yet. 

Let us now dwell upon the compatibility of the fit value for g{ J /if}) with the concept of 
J /if) — u) — (j> mixing. In this approach one would write the reduced amplitude as 

5 mix ( J/ip) ~ Mf ^ 6 C sin By g^o (s = M^) . (14) 

The Lop°7T° vertex function at s = can be estimated from chiral anomaly predictions for the 
coupling constants <7p 77 and vector meson dominance^} 

6M p M w 

g wn o n o(s = 0) ~ — - — . (15) 

Using our values for the J /if; — uj mixing parameters (p0|) and a monopole form for the s 
dependence of the vertex function, A 2 /(Mj^ — A 2 ), we find that the fit value of g mix (J/ip) 
(|T3| ) is reproduced for A ~ 1 GeV. With regard to all the uncertainties encountered in an 
estimate like this, a value of about 1 GeV for A appears reasonable and is compatible with the 
typical inverse radius of a light hadronic system. Thus, we conclude that the small probability 
of the light-quark Fock components of the J /if) as estimated by J /ip — to — <p mixing, suffices to 
generate the large J/ip — * PV branching ratios seen in experiment since it is overcompensated 
by the very large soft u,<f>—* PV transition. 



7 This approach has been used to predict the radiative PV transition, see Ql an d references therein. Kramer, 
Palmer and Pinsky investigated the u>p°7r° vertex on the basis of effective Lagrangians. Their numerical results 
are in agreement with Eq. (|l5|). They also found \g LU ' p a n o(s = 0)| <C |i?a;p 7r ( s = 0)| which is in line with our 
assumption <? mix (*') < ff mix (J/V'). 
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vv 


Cyv 






pp 


1 


16 


25 ±8 p 




LOU) 


1 


5.2 


< 3.1 g 










< 6.3 gj 




K*W 


(i/ + V2cat(&>-0 8 ))/2 


11.1 


8.2 ±3.9 [ 


1 








9.0 ±5.0 


1 


H> 


yV2cot{(l>p-9 8 ) 


1.4 


7.1 ±2.2 


g 








3.1 ± 1.2 [ 


if 



Table 5: Flavor factors (quoted for single channels) of decay amplitudes as well as experimen- 
tal and theoretical branching ratios for r/ c — * VV summed over all charge states. The quoted 
theoretical values correspond to a fit with g mix (r] c ) = 0.26 GeV; the mixing parameters for 
the pseudoscalar mesons are 4>p = 39.3°, 9$ = —21.2° and y = 0.81 pCp . 

4 Decays rj c — ► VI/ 

We write the helicity amplitudes for the r/ c — ► decays in a similar fashion as in Eq. (Q) 
and parametrize the invariant amplitude analogously to Eq. ( |l0|) as 



Ayyfe) = c^^-^), (16) 

with the decay width given by 

^-^-iSF^lW. (17) 

We assume that the r/ c decay is mediated by the mixing of the rj c meson with the r\ and 77'. 
The ij — rj' mixing parameters have already been quoted. The mixing mechanism leads to the 
flavor factors Cyy shown in Table [| (c^o 1 ^ = !)• Electromagnetic contributions to r] c decays 
can only proceed via two photons and are thus negligible (see also footnote ||) . Also a doubly 
OZI-rule violating decay mechanism should be strongly suppressed here since the production 
of light vector mesons through gluons is not enhanced by the U{1)a anomaly. The possibility 
of a decay through the cc Fock component of the light mesons is, as in the case of the J/t/i, 
ignored, see the discussion in Sec. 0. 



A fit to the experimental branching ratios fj [^, 47, 48] yields g mix (r] c ) = 0.26 GeV. The 
corresponding results are listed in Table |5[ The strength of g mix (r] c ) is compatible with the 
concept of mixing. To see this we first note that here the same vertex function, gpvv, enters as 
in the J ftp — > PV decays in the vector meson mixing approach. Analogously to Eqs. ( p^| , p^ ) 
the reduced amplitude reads 

5 mix ( % ) ~ Ml\9 c \ S m(cj>p-6 s )g Vqp o p0 (s = Ml), (18) 

and the vertex function at s = is estimated by 



W(' = Q)* wit/f • (19) 



8 Since the two results for the </><j> channel quoted in Refs. ^] are not compatible with each other, we 
have doubled the respective errors in the fit. 
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Inserting the value - 
the r\ c as well as f q 
function the value g mv 



1.0° for the mixing angle 9 C that controls the light quark admixture to 



1.07 f n [2C] and assuming again a monopole behavior of the vertex 
(Vc) = 0-26 GeV is reproduced for A ~ 1.3 GeV. This indicates a 
milder s dependence of the vertex function in Eq. ( |l~8|) than in Eq. (14). We stress that the 
mixing approach explains the order of magnitude of the rj c — > VV decay widths correctly. 
This is, to our opinion, a highly non-trivial fact. In tendency our results agree better with 
the experimental results quoted by the DM-2 collaboration than with the ones found by 
MARK-III [47, 48]. Better data for the rj c — ► VV channels are obviously needed. In analogy 
to the ip' — ► PV decays the r]' c — > VV decays are expected to be strongly suppressed. 



5 Comments on decays into baryon-antibaryon pairs 

A perturbative calculation of the J/ip and ^f' decays into baryon-antibaryon pairs has been 
carried through in Ref. j||] using the light-cone wave function for baryons proposed in 
Ref. pL9[| . In contrast to the Xcj decays here it suffices to consider the color-singlet con- 
tributions only. The color-octet contributions are suppressed by powers of l/m c and start at 
the order of the charm-quark velocity, v, squared. Moreover, there is no obvious enhancement 
of the corresponding hard scattering amplitudes which appear with at least the same power of 
a s as the valence Fock state contributions. Thus, despite the fact that m c is not very large and 
v not small (v 2 ~ 0.3), it seems reasonable to expect small higher Fock state contributions to 
both these decays. One may, however, wonder what the strength of the mixing mechanism is. 
Since this mechanism is responsible for the large PV branching ratios (see Sec. Q) it may also 
contribute to the baryon-antibaryon channels substantially and spoil the results presented 
in Ref. [42]. An estimate of the mixing contribution in analogy to Eqs. ( p^Jl5| ), using, for 



instance, the w-nucleon coupling constants 5y ajp p/4vr = 8.1 ± 1.5 and gf^pp/^ 71 " = 0-14 ± 0.2 



[50], reveals that it amounts to less than 20% of the experimental width if the s dependence 
of the ujpp vertex function is at least as strong as that of the ujpn one. Since the s dependence 
of a baryon vertex is likely to be stronger than that of a meson vertex the mixing contribu- 
tion to the baryon-antibaryon channels is probably very small. Therefore the J/ip and VP' 
decays into baryon-antibaryon pairs seem to be under control of perturbative QCD. Indeed, 
the order v° perturbative results obtained in Ref. [42] agree well with experiment [||, pi]]. As 
we emphasized in Sec. § the dominance of the perturbative contribution is further supported 
by the fact that the VP' — J/ip ratio K^-g, defined in Eq. @, is in agreement with unity (see 
Table §) as is indicative for short-distance cc annihilation. 

The hadronic helicity non-conserving decay r/ c — ► pp may be estimated through r] c — rj — rj' 



mixing. Using the coupling constant g 2 — /Air ~ 1 and g , N jr ^ [45, 50] and equating 



the s dependence of the rjNN vertex with that of the rjVV one as determined from the fit 
discussed in Sec. we find B(r] c — > pp) ~ 0.5 • 10~ 3 . Comparing with the experimental value 
of (1.2 ± 0.4) • 10~ 3 H we are tempted to conclude that the mixing mechanism also controls 
the rj c — > pp decays. 



6 Concluding remarks 

Exclusive charmonium decays constitute a laboratory for investigating power corrections and 
higher Fock state contributions as well as for studying the interplay of pQCD and soft mech- 
anisms. That there are two distinct dynamical mechanisms at work can most easily be seen 
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by comparison of J/tp — > PV and VP' — ► PV decays. Their scaling with the charmonium 
decay constants is indicative for a short-distance mechanism (as, for instance, cc annihila- 
tion into virtual gluons or photons) while a clear deviation from that scaling (as is seen for 
the helicity non-conserving PV channels) signals the prominent role of a more peripheral 
soft mechanism. Several of the charmonium decays into pairs of pseudoscalar and/or vector 
mesons are forbidden by spin and parity invariance. Others, for instance the decays of the 
C-even charmonia into PV channels, are under control of the electromagnetic decay mech- 
anism because either G-parity or isospin conservation is violated for the channels involving 
non-strange mesons. Therefore, only the following classes of decays are actually possible: 
the helicity non-conserving J/tp{^') — > PV and rj c — > VV decays which are not dominated 
by leading-twist pQCD; and the Xcj — ► PP, VV (J = 0, 2) decays which are accessible to a 
perturbative treatment and the J/ip(^f') — > PP, VV decays. The latter decays are of electro- 
magnetic nature, the amplitudes being related to the time-like PP or VV electromagnetic 
form factors at s = M 2 (n 3 Si). One should, however, be aware of corrections to the elec- 
tromagnetic decay mechanism from the isospin-violating part of pQCD. In principle, allowed 
by strong interactions are also the helicity non-conserving Xcl ~^ VV decays. We however 
expect very small branching ratios for these channels. The decays of the C-even charmonia 
into the PV channels proceed through cc annihilation into two photons and are therefore 
strongly suppressed. In fact these decays as well as the Xcl — ► VV have not been observed in 
experiment as yet ||. 

The perturbative analysis of the % c j — ► PP decays has been carried through in Refs. [||, |23|] 
and the results have been found to agree fairly well with experiment although, with regard 
to the recent BES data |3^], some readjustment of the wave function parameters seems to 
be advisable. The VV channels, on the other hand, have not been analyzed as yet. Note 
that the required light-cone wave functions of the vector mesons are not known to a sufficient 
degree of accuracy. 

In the present paper we focus on the J/ ip and VP' decays into the PV channels as well as the 
r] c — ► VV decays. We assume that in these cases the charmonium state decays dominantly by 
a soft mechanism through Fock components built up from light quarks only. This mechanism 
probes all quark-antiquark separations in the charmonium wave functions and therefore feels 
the difference between a IS and 2S-wave function. We model this mechanism by J/ip — uj — 4> 
and rj c —rj—rj' mixing, respectively, and argue that this contribution is negligible for VP' decays. 
The light vector and pseudoscalar mesons couple through vertex functions to the final state 
mesons. With a few parameters adjusted to experiment we find fair agreement with the many 
channels measured for J/?/>(\P') — * PV and rj c — * VV decays. Most of the parameters have 
a simple, direct physical interpretation and their magnitudes can be estimated from physical 
considerations. Data of improved quality, in particular for the VP' and r] c decays are required 
for a more stringent test of our approach. A richer database for the vector-tensor channels 
in J ftp and VP' decays would also allow the extension of our mixing approach to this class of 
reactions. 
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Appendix: Vector meson mixing 



In order to estimate the strength of the light-quark Fock components in the J ftp we consider 
J /ijj — u — 4> mixing in parallel to the mixing of the pseudoscalar mesons. As in Refs. [^C], 37] 
(see also f45fl ) we start from the quark-flavor basis formed by the three states u q , u> s and u c 
which are assumed to have the parton decompositions 



q/ 



^ q | UU+dd)/ V / 2 + 

* s | ss) + • • • , 

*c I CC) + • • • 



(20) 



where the ellipses stand for higher Fock states. The physical meson states are related to the 
basis ( p0|) by an orthogonal transformation 



U) 



J/4> 



U(</>v,< 



y, v C ) 



( 



(21) 



where 



U 



I 
\ 



1 

■ sin( 



-cpv 
1 

. COS I 



-# c sin l 
8 C cos i 
1 



(22) 



(Wll = 1 + 0(9%, 4>y, c 4>v))- In Eq. ( p2|) 9 C and 4>v are considered as small quantities. 9 C is 
small since the mixing between the light and the heavy sector is an effect of 0(1/ ml) while 
the smallness of <j)y follows from the well-known fact that the u and the <f> are nearly ideally 
mixed 51]. 

The mass matrix in the quark-flavor basis is related to the (diagonal) mass matrix for the 
physical mesons by 



M 



qsc 



tf(ct> v ,9 y , 9 C ) diag[M2, Ml M 2 J/f ) U^v, y , 9 C ) . (23) 

, we parametrize the mass 



Following our previous work on the pseudoscalar mesons 
matrix as 



/ 



M 



qsc 



m\ q + 2a 2 



yV2a 



7,2 



yV2a 2 
m 2 ss + y 2 a 2 
yza 2 



7.2 



yza 2 

2 

ml. 



(24) 



where y and z allow for flavor symmetry breaking effects in the otherwise "democratic" mixing 
matrix whose strength, a 2 , being controlled by quark-annihilation processes. 

Comparing the expressions (23) and ( p4|) we find the following six relations between the 
nine parameters 



m 2 qq + 2d 2 



l ss 
2 



^a 2 

r'2 



Ml , 



Mi 



1 



M 



V2ya 2 



j/i> > 
(Ml 



17 



tan 9„ 



2y , 



2 + pza 2 /M^ 



(25) 



One way to utilize these relations is to make use of flavor symmetry in order to fix the quark 
mass terms. Thus, one may use 



m\ q ~ 



M 



p ' 



2M 2 K * 



M 



p ' 



(26) 



in analogy to the pseudoscalar meson case [2C]. This would allow the determination of the 
mixing parameters with the exception of z. However, we dismiss this method since it leads 
to values for the parameters which are instable under small corrections to the relations (p6|). 

An alternative method that leads to robust values for the mixing parameters, is to take 
the mixing angle <j)y from experiment and to fix the parameters controlling flavor symmetry 
breaking in analogy to the pseudoscalar case by 



y = fu/U, z = fu>/fj/i> 



(27) 



This is an assumption here in contrast to the pseudoscalar case where the divergences of the 
axial vector currents, embodying the U{\)a anomaly, relate the flavor symmetry breaking 
parameters to the meson decay constants. From the values of the decay constants quoted in 
Ref. |H, we obtain y = 0.80 ± 0.04 and z = 0.48 ± 0.02. The mixing ansatz (H) directly 
relates the ratio of the widths for the radiative 
angle <p v || 

T(<f> - A) 



7T and U) — 7r transitions to the mixing 



T(u> — > 7T°7) 

From the experimental decay widths § we find 



Mi 

Mu. 



(28) 



<j> v = 3.4° ± 0.2 C 



(29) 

where the sign is chosen according to Eq. (]25|). This result is consistent with other determi- 
nations of that angle 52 ] . 

Using Eqs. (p7| , p9| ) and of course the masses of the physical vector mesons as input to 



Eq. (25), we obtain for the other mixing parameters 



a 2 



0.022 ±0.002 GeV 2 
-59.8° ± 1.3° , 
-0.10° ±0.01° . 



(30) 



The mixing strength, a 2 is much smaller than in the case of the pseudoscalar mesons 
where it is dominated by the U{V)a anomaly. For completeness we note that an evaluation 
of the quark mass terms from these values for the mixing parameters leads to values which 
merely deviate by about 2% from those obtained via the flavor symmetry relations (l2f 
Our approach leads to the following quark content of the physical mesons 



Wo 



0.060 \u n 



\lj s ) + 0.060 \u 



1.5 • 10~ 3 
0.9-10"" 



\u c ) + 1.5 • 10 



-3 



+ 0.9 • 10~ 3 \uj s 



(31) 
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